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Abstract 

Nowadays control systems are mostly implemented on digital plat- 
forms and, increasingly, over shared communication networks. Reducing 
resources (processor utilization, network bandwidth, etc.) in such imple- 
mentations increases the potential to run more applications on the same 
hardware. We present a self-triggered implementation of linear controllers 
that reduces the amount of controller updates necessary to retain stabil- 
ity of the closed-loop system. Furthermore, we show that the proposed 
self-triggered implementation is robust against additive disturbances and 
provide explicit guarantees of performance. The proposed technique ex- 
hibits an inherent trade-off between computation and potential savings 
on actuation. 

1 Introduction 

The majority of control systems are nowadays implemented on digital platforms 
equipped with microprocessors capable of running real-time operating systems. 
This creates the possibility of sharing the computational resources between con- 
trol and other kinds of applications thus reducing the deployment costs of com- 
plex control systems. Many control systems are also implemented over shared 
communication media making it necessary to share the communication medium. 
The concept of self-triggered control was introduced by Velasco et al in [T] to 
take advantage of the possibility (or necessity) of sharing resources. The key 
idea of self-triggered control is to compute, based on the current state measure- 
ments, the next instant of time at which the control law is to be recomputed. 
In between updates of the controller the control signal is held constant, and 
the appropriate generation of the update times guarantees the stability of the 
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closed-loop system. Under a periodic implementation, the control law is exe- 
cuted every T units of time, regardless of the current state of the plant. Hence 
this period T has to be chosen in order to guarantee stability under all possible 
operating conditions. On the other hand, under self-triggered implementations 
the time between updates is a function of the state, and thus less control ex- 
ecutions are expected. On the other hand, the intervals of time in which no 
attention is devoted to the plant pose a new concern regarding the robustness 
of self-triggered implementations. 

The contribution of this paper is to describe a self-triggered implementa- 
tion for linear systems, in which the times between controller updates are as 
large as possible so as to enforce desired levels of performance subject to the 
computational limitations of the digital platform. By increasing the available 
computational resources, the performance guarantees improve while the number 
of controller executions is reduced. Hence, the proposed technique reduces the 
actuation requirements (and communication, in networked systems) in exchange 
for computation. Furthermore, we also show that the proposed self-triggered im- 
plementation results in an exponentially input-to-state stable closed-loop system 
with corresponding gains depending on the available computational resources. 
A preliminary version of these results appeared in the conference papers [5] 
and 

Several self-triggered implementations have been proposed in the last years, 
both for linear [1] and non-linear [5] plants. The latter when applied to lin- 
ear systems degenerates into a periodic implementation, while [4] makes use of 
very conservative approximations. In contrast with those two techniques the 
approach followed in the present work provides large inter-execution times for 
linear systems by not requiring a continuous decay of the Lyapunov function 
in use, much in the spirit of [5]. Computing exactly the maximum allowable 
inter-execution times guaranteeing stability requires the solution of transcen- 
dental equations for which closed form expressions do not exist. Our proposal 
computes approximations of these maximum allowable inter-execution times 
while providing stability guarantees. The idea advocated in this paper, trading 
communication/actuation for computation, was already explored in [7j. How- 
ever, their approach is aimed at loosely coupled distributed systems, where local 
actuation takes place continuously and communication between subsystems is 
reduced by means of state estimators. In the analysis of robustness of the pro- 
posed implementation the authors were influenced by the approach followed 
in [8] and ^ . Finally, the notion of input-to-state stability [10] is fundamental 
in the approach followed in the present paper. 

2 Notation 

We denote by the positive real numbers. We also use — M+ U {0}. 
The usual Euclidean (^2) vector norm is represented by | • |. When applied to a 
matrix | • | denotes the I2 induced matrix norm. A matrix P S K™^™ is said to 
be positive definite, denoted P > 0, whenever x^Px > for all x ^ 0, x G R"*, 



2 



and a matrix A is said to be Hurwitz when all its eigenvalues have strictly 
negative real part. We denote by / the identity matrix. By Xm{P), ^m{P) we 
denote the minimum and maximum eigenvalues of P respectively. A function 
7 : [0,oo[— >■ M.Q , is of class /Ceo if it is continuous, strictly increasing, 7(0) = 
and 7(s) — > cxd as s — >■ 00. Given an essentially bounded function 5 : Rq — > M"* 
we denote by ||5||c>o its jCqc norm, i.e., ||^||oo — (ess) supjgjj+ {\S{t)\} < 00. We 
consider linear systems described by the differential equation — + Bx + S 
with inputs x • ^0 ~^ ^iid ^ • ^0 ~^ essentially bounded piecewise 
continuous functions of time. The input x will be used to denote controlled 
inputs, while d will denote disturbances. We refer to such systems as control 
system,s. Solutions of a control system with initial condition x and inputs x ^-nd 
(5, denoted by ^a:x5> satisfy: £,xx&{^) = ^ and -^^xxsit) = A£,:oxs{t) + Bx{t) + 6{t) 
for almost all t G Mg . The notation will be relaxed by dropping the subindex 
when it does not contribute to the clarity of exposition. A linear feedback law 
for a linear control system is a map u = Kx; we will sometimes refer to such a 
law as a controller for the system. 

Definition 2.1 (Lyapunov function). A function V : M™ — >■ M.q , is said to be 

a Lyapunov function for a linear system ^ = A£^ if^x^ K™ a(|a;|) < V{x) < 
a(|a:|) for some a,a€ K.00, and there exists A G R"*" such that for every x € M™; 

dV 

^Ax < -XVix). 

ox 

We will refer to A as the rate of decay of the Lyapunov function. In what 
follows we will consider functions of the form V{x) = {x^ Px) 2 , in which case V 
is a Lyapunov function for system ^ = if and only if P > and A^P + PA < 
—2X1 for some A G M+, the rate of decay. 

Definition 2.2 (EISS). A control system ^ = A^ + 6 is said to be exponentially 
input-to- state stable (EISS) if there exists X G R+, a G 1R+ and 7 e ICoo such 
that for any t e and for all x e : 

\U{t)\<(T\x\e-^'-\-j{\\S\U. 

We shall refer to {13, ■j), where /3(s,<) = sae~'^*, as the EISS gains of the 
EISS estimate. If no disturbance is present, i.e., 5 = 0, an EISS system is said 
to be globally exponentially stable (GES). 



3 A self-triggered impleraentation for stabiliz- 
ing linear controllers. 

Consider the sampled-data system: 

m = Am + Bxit) + Sit) (1) 

X{t) = K^{tk), t€[tk,tk+i[ (2) 
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where {tfc}j.gf^ is a divergent sequence of update times for the controller, and A+ 
BK is Hurwitz. The signal 5 can be used to describe measurement disturbances, 
actuation disturbances, unmodeled dynamics, or other sources of uncertainty as 
described in 1^. 

A self-triggered implementation of the linear stabilizing controller ^ for the 
plant (H)) is given by a map F : M™ — )■ R+ determining the controller update time 
tk+i as a function of the state £,{tk) at the time i^, i.e., tk+i = ifc + r(^(ifc)). 
If we denote by the inter-execution times Tk = tk+i — t^, we have = 
r(^(tfc)). Once the map F is defined, the expression closed-loop system refers 
to the sampled-data system ([Ij and © with the update times tk defined by 

tk+i = tk + r{c{tk)). 

The problem we solve in this paper is the following: 

Problem 3.1. Given a linear system {Ip and a linear stabilizing controller 
construct a self -triggered implementation F : R™ — >■ R"*" of (0) that renders EISS 
the closed-loop system defined by (QP, while enlarging the inter- execution 
times. 

In order to formally define the self-triggered implementation proposed in this 
paper, we need to introduce two maps: 

m he, a continuous-time output map and 
• /id, a discrete-time version of he. 

Let be a Lyapunov function of the form 

V{x) — (x'^Px)^ for £^ = {A-\- BK)^, with rate of decay Aq. The output map 
he : R™ X R(]" ^ R(}- is defined by: 

he{x,t):=V{Ut))-V{x)e-^' (3) 

for some < A < Aq. Note that by enforcing: 

he{Utk).t)<0, Vt e [0,T,[VfceN (4) 

the closed-loop system satisfies: 

V{£,x{t))<V{x)e-^\ VteR+VxeR" 

which implies exponential stability of the closed-loop system in the absence of 
disturbances, i.e., when 5{t) = for all t S R,^. 

Our objective is to construct a self-triggered implementation enforcing 
Since no digital implementation can check ^ for all t € [tk,tk+i[, we consider 
instead the following discrete-time version of @ based on a sampling time 
A e R+: 

hd{CAtk),n) := he{Utk),nA) < Vn G [o, [^j [, 

and for all k G N. This results in the following self-triggered implementa- 
tion, first introduced by the authors in [2], where we use N^iin '■— [''min/AJ, 
^max := [■''max/AJ, and Tniin and Tniax are design parameters. A similar ap- 
proach was followed in 11 in the context of event-triggered control. 
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Definition 3.2. The map : M" ^ M+ defined by: 



Td{x) :== max{T,„i„, rifcA} 

rifc max{n < iVmax|/id(a;, s) < 0,s = 0, . . . ,n} 

neN 

prescribes a self-triggered implementation of the linear stabilizing controller (0) 
for plant {Ip. 

Note that the role of Tmin and Tmax is to enforce exphcit lower and upper 
bounds, respectively, for the inter-execution times of the controller. The upper 
bound enforces robustness of the implementation and limits the computational 
complexity. 

Remark 3.3. Linearity of fip and 0) enables us to compute hd as a quadratic 
function of ^(ifc). Moreover, through a Veronese embedding we can imple- 
ment the self -triggered policy described in Definition \3.2\ so that its computation 
has space complexity q g^^g time complexity q+ (2(7 + I'^min^tll where 
q := iVjuax ~ -^min • For reasons of space we omit these details. They can be 
found in J^. 



4 Main results 

The proofs of all the results reported in this section can be found in the Ap- 
pendix. The following functions will be used to define ElSS-gains: 

PP '■= X tnx ' 7p,t(s) := s— / |e \dr. 



XmiP) 



\7n{P) 



We start by establishing a result explaining how the design parameter Tmin 
should be chosen. The function can be seen as a discrete-time version of the 
function : M™ ^ defined by: 

Tc{x) := max{T < T,„ax|^c(a;, s) < 0,Vs G [0, r]}. (5) 

reR+ 

If we use Fc to define an ideal self-triggered implementation, the resulting inter- 
execution times are no smaller than t^-^^ which can be computed as detailed in 
the next result. 

Lemma 4.1. The inter- execution times generated by the self-triggered imple- 
mentation in are lower bounded by: 

r^i„ = min{r e M+ : det M{t) = 0} (6) 

where: 



M{t) := C{e^ ''C^PCe^^ -C^PCe-^'')C'^, 

A+BK BK 
-A - BK -BK 



F 
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The computation of t^^^ described in Lemma HTTl can be regarded as a formal 
procedure to find a sampling period for periodic implementations (also known 
as maximum allowable time interval or M ATI) . It should be contrasted with the 
frequently used ad- hoc rules of thumb [T^], [13] (which do not provide stability 
guarantees). Moreover, an analysis similar to the one in the proof of this lemma 
can also be applied, mutatis mutandis, to other Lyapunov-based triggering con- 
ditions, like the ones appearing in [14] and [4]. Notice that the self-triggered 
approach always renders times no smaller than the periodic implementation, 
since under a periodic implementation the controller needs to be executed every 
T^ijj (in order to guarantee performance under all possible operating points). 

The second and main result establishes EISS of the proposed self-triggered 
implementation and details how the design parameters Tmin, Tmax, A, and A af- 
fect the ElSS-gains. 

Theorem 4.2. // Tmin < '''min' self -triggered implementation in Defini- 
tion \3.2\ renders the closed-loop system EISS with gains (/3,7) given by: 

(3{s,t) ■— ppg{A,Nnia.x)e'^*s, 

f^ I , A~^(P)g(A,7Vn,ax) , 

7(s) := 7P,iv„axA(s ., . ^ -H7/,iv„axA(s) 

1 — e -^^mm 



where: 



(p + 2A)f.A 



(p+2A)^A 2A^A 



-f e """^ ye i^-p — e 1^- 

p ~ Am(G), fi := A„i(G), 

piAP-^-^^{P^AP~if P^BKP- 
(piBKP-i)'^ 



G ■- 



Note that while r„iin is constrained by r^jn, Tmax can be freely chosen. How- 
ever, by enlarging Tmax (and thus Nmax) we are degrading the ElSS-gains. It is 
also worth noting that by enlarging Tmax one can allow longer inter-execution 
times, and compensate the performance loss by decreasing A, at the cost of 
performing more computations. 

Let us define the maximum exact inter-execution time from x as t*{x) := 
min{rc(a;), Tmax}, where the upper bound is required to obtain robustness against 
disturbances. The third and final result states that the proposed self-triggered 
implementation is optimal in the sense that it generates the longest possible 
inter-execution times given enough computational resources. Hence, by enlarg- 
ing the inter-execution times we are effectively trading actuation for computa- 
tion. The proof of the following proposition follows from the proof of Theo- 
rem g]!] 

Proposition 4.3. The inter- execution times provided by the self -triggered im- 
plementation in Definition \S.2\ are bounded from below as follows: 



Td{x) > T*{x) - A, Va; e 
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Note that even if Tdix) > t*{x) the performance guarantees provided in 
Theorem 14:21 stih hold. 

Remark 4.4. When implementing self -triggered policies on digital platforms 
several issues related to real-time scheduling need to be addressed. For a discus- 
sion of some of these issues we refer the readers to U5\/ . Here, we describe the 
minimal computational requirements for the proposed self -triggered implemen- 
tation under the absence of other tasks. There are three main sources of de- 
lays: measurement, computation, and actuation. Since the computation delays 
dominate the measurement and actuation delays, we focus on the former. The 
computation of Td is divided in two steps: a preprocessing step performed once 
by execution, and a running step performed n times when computing hd{x,n). 
The preprocessing step has time complexity {m'^ +m)/2 and the running step 
has time complexity m? -\- m. If we denote by Tc the time it takes to execute an 
instruction in a given digital platform, the self -triggered implementation can be 
executed if: 

-{m^ + m)Tc < Tniin, (m^ + m)Tc < A. 

The first inequality imposes a minimum processing speed for the digital platform 
while the second equality establishes a lower bound for the choice of A. 

We refer the interested reader to [2] and ^ for numerical examples illustrat- 
ing the proposed technique and the guarantees it provides. 

5 Conclusions 

This paper described a self-triggered implementation of stabilizing feedback con- 
trol laws for linear systems. The proposed technique guarantees exponential 
input-to-state stability of the closed-loop system with respect to additive dis- 
turbances. Furthermore, the proposed self-triggered implementation allows the 
tuning of the resulting performance and complexity through the selection of 
the parameters A, A, r^-m and Tmax- The performance guarantees can be im- 
proved and the inter-execution times enlarged by increasing the computational 
complexity of the implementation. 
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6 Appendix: Proofs 

Proof of Lemma\4-1\ It can be verified that he satisfies hc{x^{i) ~ and Mr hc{x,t) < 

" 1=0 

0,Vx e M™, which, by continuity of /ic, implies the existence of some t*^ii^{x) > 
such that Tc{x) > r^i„(a;). Let us define the variables 7]{t) = ^{t) — £,{tk), 
t € [tk,tk+i[ and C = [C"^ V^]'^ ■ Note that at the controller update times = 
0. Under this new notation, system ([T]) with controller ([2|), in the absence of 
disturbances, can be rewritten as ({t) — F({t) with solution Cy{t) = e^*y, 
y = [a;T 0^]"^. Let us denote by the map hdy^t) = V{CCyit)) - V{Cy)e-^*. 

While it is not possible to find Fc in closed form, we can find its mini- 
mum value by means of the Implicit Function Theorem. Differentiating (f>{x) = 
hc{y, Fc(a;)) — with respect to the initial condition x we obtain: 



dTc dhc dy 

— - H -— = 0. 

t=Ta(x) dx dy dx 



dcf) dhc 
dx dt 

The extrema of the map Fc are defined by the following equation: 

-1 

/ r)/? . dm \ 

0. 



dx 



dhc dy 

t=Vc{x) ) \ dy dx 



dha 
dt 



Hence, the extrema of Fc satisfy either ^ s ^ ^ for some t e M+ or ^ (x) = A 
— 0. The latter case corresponds to situations in which for some x 

t=T^(x) 

the map he reaches zero exactly at an extremum, and thus can be disregarded 

dhc dy / 
dy dx ^ 



as violations of the condition hc{tk,t) < 0. Combining ^^^{t,x) = into 
matrix form we obtain: 

M{t)x = 0. (7) 

The solution to this equation provides all extrema of the map Fc(a;) that incur 
a violation of hc{x,t) < 0. Thus, the minimum r satisfying 1^ corresponds to 

d 



the smallest time at which hc{x,T) — 0, -^hdx^t) 



> can occur. Since the 

left hand side of ([7|) is linear in x, it is sufficient to check when the matrix has 
a nontrivial nullspace. Hence the equality ([6]). □ 



We introduce now a Lemma that will be used in the proof of Theorem 14. 2 1 
Lemma 6.1. Consider system (Qp and a positive definite function V{x) = 
(x'^Px)^ , P>Q. For any given < T < oo the following bound holds: 

v{^xx5{t)) <vi^^^oit)) + -fP,Tm\o.),yt e [o,T]. 

Proof. Applying the triangular inequality and using Lipschitz continuity of V 
we have: 

XUP) 
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Integrating the dynamics of ^ and after applying Holder's inequality one can 
conclude that: 



t 

Ar 



%At) - i.xo{t)\ < / |e^1rfr||<5||oo. 
And thus for aU t e [0,T]: 

A7„(P) Jo 



□ 



Proof of Theorem \4-S\ We start by proving that in the absence of disturbances 
the following bound holds: 

\Utk+r)\<g{A,Nn,.^)\Utk)\e-^\ Vr > 0. (8) 

Let W{x) — x^Px and use W{t) to denote W^(Cxfc(i)), with determined 
by (P), ([2]), and = ^di£.{tk))- By explicitly computing W{t) one obtains: 

\pk{t)f {piah)r]G[{piat)r {pH{tk)fY , 

for t £ [tk,tk+i[, and thus the following bounds hold: 

M {Wit) + W{tk)) < W{t) < p {W{t) + W{tk)) . 

for t £ [tk,tk+i[- After integration, one can bound the trajectories of W{t), 
when t + s belongs to the interval [tk,tk+i[, as: 

W{t + s) < e"''W{t) + W(tk){eP'' -1), 
W{t + s) > e'''W(t) + W{tk){e'"' -1). 

Let us denote tk + nA by r„ for succinctness of the expressions that follow. An 
upper bound fW(W^(+)s^^d for t e [r„,r„+i[ is then provided by: 

( e'"'Wir„) + W{tk){e>" -1), se[0,s'] 

i e*^"^""^' W"(r„ + A) + W'(tfc)(e'''*"^-' — 1) S £ [s* Af. 
The maximum for the bound^of (r„ + s) for s G [0, is Attained at the 

point at which the two branches of the bound meet, i.e. at s = s*, as the first 

branch is monotonically increasing in s, and the second branch monotonically 

decreasing. The point s* can be computed as: 



p-fi ""V W{rn) + W{tk) 
and thus W{rnW($*<^^m)\^ bSKi^M^dl-as: 
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which is monotonically increasing on W(r„), VK(r„+i), and W{tk). If S{t) = 
W(tfc)e-2^(*-*'=),H{(fcatews) < -S(tfc) + 

where we used the raci that, ii Tmin < '''mim ^ d enforces (in the absence 
of disturbances) W{rn) < S'(r„) for all n G N, n < n^. From the previous 
expression we can obtain W(r„ + s*) < g{A,n)S{rn + s*) where: 



f/(A, n) — — e ^ + 
e'^-Ple (1 + e )l 



jrr^/ -2A(A-s*) _|_ ^2A(„A + s')N JT^ 



The value of s* can be further bounded to obtain a simpler expression: 

/i - p 

Using this bound for s* and letting n take its maximum possible value n = 
-^max — 1, the following chain of inequalities holds: 

pp.9(A,n)3 < pp5(A,7V,nax- 1)^ < <7(A,iVmax) 

for all n G [OjiVmax], which leads to the bound: 

W^-{t) < pp'giA,N^,^)sHt). (9) 
Note that ([9]) does not depend on or n. Finally, apply the bounds: 



XliP)\x\ < Vix) = Vx'^Px < XIj{P)\x\. (10) 

to obtain From Lemma I5TT1 and the condition enforced by the self-triggered 
implementation we have: 

Iterating the previous equation it follows: 

fe-i 



Vmk)) < e-^(*-*°)nC(0)+7p(ll'5||oo)Ee 

oo 

< e-^(*-*<')F(e(0)+7p(ll'5||oo)E 



oo 

4=0 



= e-^'*'"-*°V(e(0)+7p(ll'5||oo)^ 

1 — e -^^mm 

Assuming, without loss of generality, that to ~ 0, the following bound also holds: 
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where we used (dU]). From and Lemma ISTT] one obtains: 

\U-tk + r)\ < 9iA,N,^,^)\Utk)\e-^^ + 7/(ll^l|c 
for all r e [0, A^maxA]. Combining ([TT|) and (H^]) results in: 

\Utk+T)\ < g(A,iV,„ax)pp|x|e-^(*^+^) 

, Arn^(P)g(A,iV„ 
+ 7P(II'5||«.) ^ _ g,At„„„ 

+ 7/(ll'5||oo), 

and after denoting tk + t hy t we can further bound: 
1^,(01 < g{A,N,^,^)pp\x\e-^' 



7p(ll'5|loc) ^ „ e-^*""» +7/(l|o|loo), 



which is independent of k and concludes the proof. 
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